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1  Introduction 

The  research  performed  by  our  group  at  the  Istituto  di  Struttura  della  Materia  of  the  Univer¬ 
sity  di  Messina,  Italy,  under  the  terms  of  Contract  DAJA45-86-C0003  can  be  summarized 
as  follows: 

1.  Comparison  of  the  experimental  data  on  the  scattering  properties  of  aggregated  spheres 

with  the  results  of  our  calculations; 

2.  Study  of  the  effect  of  clustering  and  correlation  on  the  optical  properties  of  intermedi¬ 

ate  and  high-density  dispersions  of  spherical  panicles  within  the  framework  of  an 
enhanced  version  of  the  Bruggeman  model. 

3.  Density  dependence  of  the  absorption  coefficient,  in  the  low  and  intermediate-density 

regime,  for  a  dispersion  of  metal  particles. 

As  a  result  of  the  research  mentioned  above  a  number  of  communications  have  been  pre¬ 
sented  to  International  Conferences,  one  paper  containing  final  results  has  been  submitted 
to  J.  of  the  Optical  Soc.  Am.  and  a  thesis  for  the  graduation  in  physics  has  been  submitted 
to  the  Faculty  of  Sciences  of  the  University  of  Messina. 

2  Comparison  with  experimental  data 

The  programs  we  developed  in  the  last  ten  years  to  calculate  the  properties  of  clustered 
spheres  proved  to  be  suitable  to  deal  with  a  number  of  physically  significant  situations. 
Nevertheless,  since  we  were  asked  in  several  instances  how  our  results  compare  with  the 
experimental  data,  we  determined  to  asses  the  reliability  of  our  calculations  through  direct 
comparison  with  the  data  for  single  clusters  published  by  Schuermann  and  Wang.  These 
authors  in  fact  performed  a  series  of  measurements  of  the  forward-scattering  amplitudes 
of  single  clusters  composed  of  2,  4  and  8  spheres  both  with  real  and  with  complex  refrac¬ 
tive  index  as  a  function  either  of  the  angle  of  incidence  of  the  incoming  plane  wave  or 
of  the  separation  of  the  components.  The  measurements  were  actually  performed  in  the 
microwave  range  on  dielectric  spheres  with  size  parameter  ranging  from  3  to  8. 

Thanks  to  the  flexibility  of  our  programs  we  were  able  to  reproduce  the  measurements 
within  the  experimental  error,  that  according  to  the  statement  of  the  authors  is  of  the  order 
of  10%  in  magnitude  and  12°  in  phase.  Furthermore,  thank  to  our  analysis  of  the  transfor¬ 
mation  properties  of  the  forward  -scattering  t'.mplitude  under  changes  of  the  direction  of 
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in  incidence  we  are  able  to  justify,  on  theoretical  grounds,  the  shape  of  tbe  curves  reported 
.  by  Schuermann  and  Wang. 

The  preliminary  results  of  our  calculations  formed  the  subject  of  a  Communication 
to  be  published  in  the  Proceedings  of  the  European  Aerosol  Conference,  held  in  Vienna, 
Austria,  on  September  18-23, 1989. 

3  Influence  of  clustering  and  correlation 

The  optical  and  dielectric  behavior  of  high-density  dispersions  of  scanerers  is  commonly 
dealt  with  within  the  frame  work  of  the  Effective  Medium  Theory  using  either  the  Brugge- 
man  or  the -Maxwell-Gamett  model.  Unfortunately,  both  these  models  were  developed 
in  the  dipole  approximation,  i.  e.  neglecting  the  effect  of  the  higher  multipoles  and  in 
particular  of  the  magnetic  dipole.  Although  in  recent  times  the  Effective  Medium  Theory 
was  modified  to  deal  with  phenomena  definitely  attributable  to  magnetic -dipole  effects, 
even  the  most  recent  improvements  are  not  applicable  but  to  particles  with  very  small 
size-parameter. 

To  improve  on  this  point  and  to  overcome  the  limitations  of  the  theory  due  to  the  com¬ 
plete  lack  of  interparticle  correlation  in  the  distribution  we  propose  some  modifications 
to  the  original  mixing  rule  of  Bruggeman.  The  first  improvement  is  to  use  the  full  Mie 
expansion  to  describe  the  the  scattering  properties  of  single  panicles:  this  allows  us  to 
deal  with  scatterers  with  size-parameters  up  to  1. 

As  a  second  improvement  we  introduced  into  a  distribution  of  dielectric  and  metal 
particles  a  correlation  of  exclusion  designed  to  prevent  contact  among  the  metallic  one. 
This  effect  is  achieved  by  putting  around  each  particle  a  thin  dielectric  layer  with  e  =  1 . 
Then,  by  gradually  increasing  the  thickness  of  the  exclusion  layer  we  were  able  to  follow 
the  transition  from  a  behavior  typical  of  the  Bruggeman  model  to  one  more  appropriate  to 
the  Maxwell-Gamett  model. 

As  a  final  improvement  we  included  into  our  calculations  the  possibility  that  the  par¬ 
ticles  aggregate  in  pairs  thus  forming  anisotropic  scanerers  either  all  oriented  alike  or  at 
random. 

The  results  we  got  can  be  summarized  as  follows.  The  effective  Medium  Theory  can¬ 
not  be  appiied  to  media  composed  of  panicles  with  size-parameters  of  the  order  unity:  our 
calculation  show  that  a  size-parameter  of  the  order  0.1  is  a  practical  limit  for  the  reliabil¬ 
ity  of  the  theory.  Furthermore,  the  clustering  may  have  visible  effects  on  the  percolation 
threshold  of  a  dispersion  of  metal  spheres,  specially  when  the  clusters  are  allowed  to  orient 
all  alike.  Finally  the  introduction  of  the  correlation  of  exclusion  produces  a  rather  sharp 
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transition  from  the  Bruggeman  model  to  the  Maxwell-Gamett  model.  Nevertheless  the 
transition  is  not  so  sharp  as  to  prevent  us  from  following  the  change  of  the  behavior  of  the 
dispersion. 

The  calculations  summarized  above  form  the  subject  of  a  Communication  to  the  2d 
Int.  Conf.  on  Optical  Particles  Sizing  held  in  Phoenix,  Arizona,  4-7  march,  1990.  A 
complete  paper  including  all  our  results  has  been  submitted  to  J.  Optical  Soc.  Am. 

4  Density  dependence  of  the  absorption  coefficient 

The  problem  of  the  density  dependence  of  the  absorption  coefficient  of  a  dispersion  of 
particles  is  of  paramount  importance  both  from  the  theoretical  and  from  the  experimental 
point  of  view.  It  is  well  known,  indeed,  that  one  often  speaks  of  low-density  or  of  high- 
density  dispersions,  but  on  the  nrecise  definitions  of  these  terms  there  exist  only  some  rule 
of  thumb  based  on  the  knowledge  gained  from  experiments.  For  this  reason  we  put  Mr 
Fucile,  one  of  our  undergraduate  students,  in  charge  of  a  research  aimed  at  developing 
a  suitable  method  to  determine  the  density  dependence  of  then  absorption  coefficient  of 
a  dispersion  of  spherical  particles,  as  pan  of  a  graduation  thesis  to  be  submitted  to  the 
Faculty  of  Sciences  of  our  University. 

The  method  used  by  Mr  Fucile  is  based  on  the  fact  that  the  density  dependence  of 
the  optical  properties  of  a  dispersion  is  due  to  the  progressive  increase  of  the  multiple 
scattering  effects  and  that  effects  of  this  kind  are  fully  accounted  for  by  the  method  we 
used  till  now  to  deal  with  clusters  of  spheres.  Therefore  the  dispersion  was  considered 
as  composed  of  pairs  of  spherical  scatterers.  The  properties  of  each  pair  were  calculated 
through  our  usual  method  including  the  analytic  average  over  the  orientations.  The  product 
of  this  calculation  is  the  average  forward-scattering  amplitude  of  a  pair  which  still  depends, 
however,  on  the  separation  of  the  component  spheres.  At  this  stage  we  take  advantage  of 
the  results  of  statistical  physics  to  set  up  the  pair  correlation  function  appropriate  for  the 
density  under  consideration  and  to  use  it  to  perform  an  average  over  the  separation  of  the 
particles.  As  a  result  we  get  an  effective  spherical  scatterer  dressed  by  the  effects  due 
to  the  presence,  in  the  dispersion,  of  other  particles.  At  this  stage  the  calculation  of  the 
optical  properties  of  the  dispersion  is  only  a  matter  of  using  well  known  formulas. 

As  noted  above,  this  research  is  the  subject  of  a  graduation  thesis  (in  italian,  of  course), 
but  we  also  extracted  a  preliminary  Communication  to  the  3d  Conference  of  the  European 
Colloid  and  Interfaces  Society  to  be  held  in  Copanello,  Italy,  next  September.  We  plan  to 
submit  a  complete  paper  after  the  graduation  of  Mr  Fucile. 
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1.  INTRODUCTION 

It  is  well  known  that  the  propagation  of  the  electromagnetic  waves  trough  a  composite  can  be  studied  by  the  effective  medium 
theory  provided  that  the  particles  be  small  enough  with  respect  to  the  wavelength  that  a  dielectric  description  of  the  composite 
as  a  whole  is  still  meaningful1  -1 .  A  possible  way  to  deal  with  this  problem  is,  indeed,  to  think  of  the  the  macroscopic  properties 
of  a  composite  as  a  space  average  of  its  local  properties  and  consequently  to  substitute  the  actual  medium  with  an  effective  one. 
In  practice  the  effective  medium  theory  achieves  this  goal  through  the  definition  of  a  homogeneous  medium  capable  of  fitting 
the  average  propagation  properties  of  the  composite  in  such  a  way  that,  if  this  latter  is  embedded  into  the  effective  medium,  no 
scattering  is  experienced1  •* .  Unfortunately  there  is  no  universal  recipe  to  build  the  dielectric  constant  of  the  effective  medium, 
since  the  properties  of  a  composite,  depend  on  the  nature,  the  shape  and  the  size  as  well  as  on  the  concentration  and  on  the 
possible  con-elation  in  the  disiribuuon  of  the  constituent  panicles. 

A  particularly  simple  model  composite  medium,  that  in  the  framework  of  the  effective  medium  theory  leads  to  the  so  called 
Bruggcman  mixing  rule,  is  composed  of  sphcncal  particles  of  very  small  size  parameter  with  totally  uncorrclatcd  space  distribu¬ 
tion.  If  into  the  model  above  is  forced  enough  corrclauon  as  to  ensure  that  the  partjclcs  arc  well  separated  and  dispersed  within 
a  homogeneous  matrix,  one  is  lead,  instead,  to  the  Maxwell-Gamed  mixing  rule2 .  Of  course,  apart  from  such  a  rough  treatment 
of  the  corrclauon.  the  approximauons  of  the  two  models  above  arc  rather  severe  as  in  actual  composite  media  the  size  parameter 
of  the  particles  is  often  not  very  small  and  the  particles  themselves  may  have  a  definitely  nonsphcricai  shape.  In  this  respect 
the  possibility  of  aggregation  of  several  spherical  particles  should  also  be  considered  as  this  phenomenon  produces  cffccuvcly 
nonsphcricai  scattcrcrs3 .  Even  by  assuming  that  all  the  particles  arc  actually  spherical,  their  radii  arc  seldom  identical  but  arc  in 
general  distributed  within  a  more  or  less  wide  range.  The  literature  reports  several  attempts  to  overcome  at  least  some  of  these 
approximations:  for  instance,  a  dynamical  extension  of  the  Bruggeman  rule  applicable  to  size  parameters  up  to  0 .01  has  been 
proposed  by  Stroud  and  Pan  1 ,  while  a  similar  extension  of  the  Maxwell -Garnett  rule  and  a  distribufion  funcuon  for  the  radii 
have  been  introduced  by  Chylck  and  Srivastava  * . 

In  this  paper  we  introduce  a  modified  Bruggeman  model  to  be  applied  when  the  size  parameter  of  the  particles  is  extended 
up  to  0  1 .  Our  extension  is  effected  by  describing  the  particles  through  the  exact  Mic  amplitudes3  •*  and  including  up  to  the 
quadrupole  terms  for  the  largest  spheres.  We  tcntauvely  considered  also  a  medium  composed  of  particles  with' size  parameter 
up  to  ~  1 .0  by  including  terms  up  to  the  2  *  -poles  to  ensure  the  convcrgcncy  of  our  calculations.  Nevertheless,  we  do  not 
regard  the  results  we  obtained  for  this  Latter  medium  as  reliable  although  they  come  out  as  a  natural  and  smooth  evolution  of 
the  results  for  spheres  with  smaller  size.  Indeed,  we  think  that  any  attempt  to  describe  the  properties  of  a  medium  too  coarsely 
grained  with  respect  to  the  wavelength  of  the  propagating  radiation  by  the  effective  medium  theory  is  bound  to  fail. 

We  apply  our  modified  model  to  study  also  how  the  correlation  affects  the  behavior  of  the  percolation  of  a  composite  made 
of  metal  and  of  dielectric  particles.  The  kind  of  correlation  we  forced  onto  this  composite  is  a  correlation  of  exclusion  among 
the  metal  particles  that  we  introduced  within  the  framework  of  our  Bruggeman-Mic  model  by  putting  around  each  metal  sphere 
a  very  ihin  dielectric  layer.  Then,  by  increasing  the  thickness  of  the  layer  we  are  able  to  follow  by  our  calculations  the  rather 
sharp  decrease  of  die  absorption  that  comes  with  the  implied  transformation  of  the  topology  from  that  of  Bruggeman  model  to 
that  of  the  Maxwell-Oamcit  one 


2.  THEORY. 

The  effective  medium  theory  defines  as  the  dielectric  constant  of  the  homogeneous  medium  that  can  be  substituted  to  the 
real  medium  without  affecting  the  propagauon  of  the  electromagnetic  waves.  This  is  accomplished  by  requiring  that  satisfies 


4 


J 


the  null  scattering  condition  and  is  thus  the  solution  to  the  equation 
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where  e,  is1  the  polarization  vector  of  the  plane  wave  and  f,ti  is  the  normalized  forward-scattering  amplitude  of  any  of  the 
identical  particles  of  the  i-th  kind  when  immersed  into  the  homogeneous  effective  medium  of  dielectric  constant  t^s- 

According  to  Mie,  the  normalized  forward-scattering  amplitude  of  a  homogeneous  spherical  particle  with  refractive  index 
n=*  y/t  embedded  in  a  homogeneous  medium  of  refractive  index  n,xl  —  in  the  present  case  r vm  =  >Vff  =  \/idi  —  is  given  by 
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where  the  polarization  index  tj  has  been  dropped  for  the  sake  of  simplicity,  k  =*  ui/c  and 
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In  cqs.  (2)  and  (3)  the  ’s  arc  the  multipolc  amplitudes  of  the  incident  plane  wave,  the  ’s  are  those  of  the  scattered  field 
and  • 
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x  a  Jcp  is  the  size  parameter  of  the  particle  with  radius  p  and  p  is  a  parity  index  that  distinguishes  the  magnetic  mullipolcs 
(p  =*  1)  from  the  electric  ones  (p  »  2).  Accordingly,  refers  to  the  magnetic  and  Rg  to  the  elcctnc  2,-polc.  When 
,|yal|  <  1  cq.  (4)  can  be  expanded  in  terms  of  ye„  and  one  gets 
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while  the  H^p) 's  for  £  >  1  are  smaller.  As  a  consequence,  in  the  limit  z*  -•  0 ,  only  the  fl*,25  's  as  given  in  cq.  (5)  can  be 
considered  in  cq.  (1)  and  one  obtains  the  well  known  Bruggcman  mixing  rule4 
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where  V,  *  A/>,/V  is  the  fraction  of  the  toial  volume  of  the  composite  filled  by  the  spheres  each  of  volume  v,,  dielectric 
constant  e*  and  size  parameter  z*.  In  order  to  get  a  dynamical  extension  of  the  Bruggcman  mixing  rule,  Stroud  and  Pan1  do  not 
neglect  the  R.\[)  term  and  thus  include  into  the  expansion  of  eq.  (2)  also  the  magnetic  dipole  term.  Since  they  take  R\l)  and 
R(,2)  in  the  form  of  cq.  (5),  their  mixing  rule  reads 
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and  thus  appears  as  an  extension  of  eq.  (6)  to  rather  larger,  although  still  small,  sizes.  The  mixing  rule  of  Chyiek  and  Srivastava4 
lies  on  the  very  same  approximations  but  these  authors  include  size-dependent  terms  only  when  e  is  very  large  —  for  instance, 
this  is  the  case  of  metal  spheres.  They  also  include  a  size-distnbution  function  X(p)  and  write 
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Obviously,  the  size-dependent  terms  in  eq.  (8)  assume  the  same  form  they  have  in  eq.  (7)  when  >,  (p)  =>  N,B{pj  -  p). 

Since  the  only  common  implementations  of  eq.  (1)  are  the  zero-size  approximation,  eq.  (6),  and  the  small-size  approxi¬ 
mation,  eqs.  (7)  and  (8),  we  determined  to  test  the  reliability  of  the  null  scattering  condition,  eq.  (1),  when  applied  to  larger 
particles.  In  practice,  all  the  calcinations  we  will  present  in  the  next  section  were  performed  by  carrying  to  full  convergency  the 
exact  Mie  expansions  to  the  effect  of  ensuring  the  inclusion  of  all  significant  contributions  to  the  scattering  amplitude,  eq.  (2). 
In  this  way  we  overcame  ail  the  limitations  due  to  the  small  size  of  the  particles  but  those  that  prevent  the  possibility  of  dealing 
with  the  whole  medium  through  a  dielectric  description. 

As  mentioned  in  the  preceding  section,  we  also  considered  a  model  composite  medium  which,  in  a  sense,  is  a  hybrid  be 
tween  that  of  Bruggeman  and  that  of  Maxwell-Gameo;  namely  we  calculated  (d;  for  a  dispersion  of  homogeneous  dielectric 
spheres  with  t  *  eD  and  of  spheres  with  a  metal  core,  of  dielectric  constant  c*,  coated  by  a  thin  dielectric  layer.  The  scattering 
properties  of  the  coated  spheres  were  not  described  through  an  approach  based  on  the  boundary  conditions  at  the  separation 
surface  but  rather  by  interposing  between  the  core  and  the  coaling  a  very  thin  transition  layer7  through  which 

<  *  «(r)  *  ty  +  (3j2  -  2 aJ)A 

with 

t  —  pm 
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where  pM  is  the  radius  of  the  metal  core  ana  po  the  inner  radius  of  the  dielectric  coating.  In  this  way  we  achieved  j  smooth 
transition  between  the  two  materials  as  the  refractive  index  n  ■  n(  r)  is  a  regular  function  of  r  everywhere  within  the  coated 
spheres.  Consequently  cq.  (4)  must  be  substituted  by  16 

where  -  G‘*\x)  arc  now  calculated  by  numerical  integration  of  die  radial  cquauons 
d7^»  2dn<  r, 

with  {  a  hr.  Of  course,  when  the  spheres  arc  homogeneous.  <o  *  and  cq.  (9)  reduces  to  cq.  (4).  As  a  matter  of  fact, 
all  the  calculations  we  will  present  in  the  next  section  were  perfotmed  by  using  cq.  (9).  or  cq.  (4)  when  appropriate,  with 
no  approximation.  It  is  also  worth  noticing  that  our  previous  calculations7  proved  that  the  introduction  of  the  transition  laver. 
provided  it  is  thin  enough,  has  no  influence  on  the  results. 


3.  RESULTS  AND  DISCUSSION. 


All  the  composite  media  we  dealt  with  in  our  calculations  arc  two-componcm  ones;  the  first  component  is  a  dielectric  with 
to  =*  1  and  the  second  is  a  metal  with 
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where  v  3  u/wp,  7  3  l/(rwp),  viz.  the  dielectric  material  is  the  vacuum  while  the  metallic  component  is  described  by  the 
frce-clectron  Drude  constant9 .  This  choice  allows  us  an  easy  comparison  with  the  preceding  results  of  Stroud  and  Pan1  and  of 
Chylck  and  Srivastava1 .  Anyway,  we  think  that  a  greater  sophistication  of  the  description  of  the  metallic  component  is  of  no 
usefulness  within  a  framework  substantially  as  naive  as  the  effective  medium  one. 

If  we  assume  the  particles  to  be  spherical,  cq.  (1)  becomes 


Vw/v(e«ff)  + (pw/po)3(l  -  Vu)fD(t«)  =0,  (11) 

where  pM  is  the  radius  of  the  metal  spheres  and  p0  th«  radius  of  the  dielectric  ones,  and  the  /’ s  are  given  by  eqs.  (2),  (3)  and 

(4). 

When  the  metal  spheres  are  coated  by  a  dielectric  layer  we  have 
VufucU* t)  +  (p«/Po)3(  1  -  Vu(l  +  cL)] f0(UH)  =*0,  (12) 


where  fuc  refers  to  a  sphere  with  a  metal  core  of  radius  p M  and  a  dielectric  coating  layer  with  external  radius 


and  is  thus  given  by  eqs.  (2),  (3)  and  (9).  Tlte  quantity  cL  is  defined  as 
cl  *  n/vu, 

where  ut  is  the  volume  of  the  dielectric  layer  in  other  words  cl  gives  the  volume  of  the  dielectric  coating  around  each  metal 
sphere  as  a  fraction  of  the  volume,  v« ,  of  the  metal  spheres  themselves.  Of  course 

^  Vurrux  *  7“ - 1  0  Cl  ^  Cl,  nus  *  Tl  1, 

1  ♦  C  L  Vy 

so  that  when  c'l  *  0 ,  i.e.  when  the  spheres  have  no  coating,  we  get  again  the  topology  of  the  Bmggernan  model,  whereas,  when 
cl  *  cL mu,  we  get  the  Maxwell-Gameit  topology. 

We  calculated  tat  according  to  eqs.  (11)  or  (12)  as  a  function  of  Vy  for  some  values  of  xy  *  kpy  and  of  u;  7  *  0.01  in  all 
cases.  In  particular  we  report  in  figs.  1  and  2  Im(  tea]  as  calculated  from  cq.  (11),  and  from  the  two-component  implementations 
of  eq.  (6),  of  eq.  (7)  with  po  *  py,  and  of  eq.  (8)  with  a  6-like  distribution  of  the  radii.  Our  main  purpose  is  to  show  the  most 
important  features  of  the  absorption  of  a  composite  medium  at  frequencies  for  which  the  metal  component  is  by  itself  strongly 
absorbing.  Therefore,  we  chose  to  plot  our  results  for  y  *  0.01,  Fig.  1,  and  for  v  ■*  0.1,  Fig.  2,  since,  for  many  metals,  these 
frequencies  occur  in  the  far  and  in  the  near  infrared,  respectively.  It  is  worth  noticing  that  both  our  eqs.  (11)  and  (12)  and  their 
counterpart  according  to  Stroud  and  Pan,  eq.  (7),  show  a  dependence  on  the  size  of  the  dielectric  spheres.  When  one  assumes 
to  *  1 ,  i.e.  when  the  dielectric  is  the  vacuum,  the  problem  arises  whether  its  more  or  less  minute  subdivision  into  spherical 
cells  may  affect  the  results.  As  a  matter  of  fact,  we  took  xp  3  kpp  *  10 -4  as  our  calculations  evidence  that  the  size  of  the 
vacuum  spheres  has  very  small  influence  cither  when  their  size  is  reduced  to  zero  according  to  Chylck  and  Srivastava  or  when 
it  is  set  equal  to  the  size  of  the  metal  spheres  according  to  Stroud  and  Pan.  Indeed,  the  results  are  nearly  indistinguishable  rpr 
all  the  sizes  of  the  metal  spheres  at  which  the  implementations  ot  the  effective  medium  theory  we  deal  with  can  be  assumed  to 
work  correctly.  The  range  of  Vy  both  in  Fig.  1  and  in  Fig;  2  has  been  chosen  so  as  to  show  the  influence  of  the  size  and  of  the 
frequency  on  the  transition  from  a  dielectric  behavior  to  a  metallic  one.  In  fact,  we  think  that,iw  =>0,1  is  an  appropriate  upper 
limit  for  the  applicability  of  cq.  (11),  viz.  it  is  a  practical  limit  to  the  size  of  the  grains  of  a  composite  medium  in  order  that  this 
latter,  from  the  point  of  view  of  the  dielectric  properties,  can  be  substituted  by  a  homogeneous  medium.  Anyway,  on  the  basis 
of  our  calculations  we  conclude  that,  for  particles  of  size  greater  that  0.1  but  rather  smaller  than  1.  the  form  of  Im(  e^J  from 
cq.  (I  t)  is  suit  meaningful  both  in  the  small- Vy  region  up  to  where  the  transition  to  the  metallic  behavior  starts  to  appear  and 
in  the  high-V*/  region.  In  this  connection  we  recall  that  cq.  (6)  should  never  be  applied  when  x  >  0 .91  and  that  eqs.  (7)  and 
(8)  cannot  be  thought  of  as  accurate  when  xy  ~  0. 1 ,  because  these  equations  come  out  from  a  senes  expansion,  cq.(5),  of  the 
exact  Mic  quantities,  cq.  (4).  We  notice  that  in  Fig.  1  the  curve  for  xy  *  0  .01  is  generated  not  only  by  our  Bruggcman-Mic 
implcmcniuuon,  cq.  ( 1 1 ),  but  also  by  the  implementations  of  Stroud  and  Pan,  cq.  (7),  and  of  Chylck  and  Srivastava,  cq.  (8);  the 
curves  from  eqs.  (7)  and  (8)  for  xy  =»  0  1  do  not  appear  in  Fig.  1  foe  the  approximations  o,.  which  eqs.  (7)  and  (8)  are  based 
result  in  solutions  that  arc  physically  unacceptable.  In  fact,  these  equations,  that  are  algebraic  equations  of  the  third  degree, 
yield  value:;  of  Imfr,*]  that  arc  cither  negauve  or  everywhere  vanishingly  small  or  very  large  even  in  the  low  Vy  region,  lh 
turn,  in  Fig  2  ihe  curves  from  eqs.  (7)  and  (8)  for  xy  *0.1  coincide  with  each  other  and  are  substantially  different  from  that 
coming  from  cq.  (11).  Moreover,  into  the  curve  labelled  (1),  which  refers  to  the  results  from  our  Bruggeman-Mie  model  for 
xy  =  0  01  .  coalesce  the  results  from  eqs.  (7)  and  (8)  for  the  same  size  as  well  as  the  results  from  the  Bruggcman  rule. 

In  Figs.  3  and  4  we  plot  Im(  t*n\  for  u  *  0  01  and  for  several  values  of  cj,  as  calculated  according  to  eq.  (14);  analogous 
results  were  obtained  for  u  *  0.1.  It  is  quite  evident  that  the  correlation  of  exclusion  produced  by  the  dielectric  coating  around 
the  metal  spheres  has  very  strong  effects  even  at  very  low  values  of  cl.  Actually,  the  model  described  in  eq.  (12)  prevents  a-.y 
coniact  among  the  metal  particles  and  is  therefore  topologically  akin  to  the  Maxwell-Garncti  model.  As  a  consequence,  it  is 
not  surprising  that,  when  cL  increases,  our  results  tend  very  fast  to  those  that  are  typical  of  the  latter  model.  Figs.  3  and  4  show 
the  results  for  ip  *  10  ~4.  At  this  value  of  xq  the  contribution  of  the  vacuum  spheres  to  the  scaucring  amplitude  is  very  well 
approximated  by  the  first  of  of  eqs.  (5)  and  such  a  small  size  suggests  the  idea  of  3  dielectric  medium  which  pervades  and  fills 
ail  Jie  space  around  the  coated  spheres.  Actually,  any  size  of  the  dielectric  spheres  up  to  xp  *  proved  to  be  incapable  of 
affecting  the  results  of  our  calculations  from  eq.  (12).  Of  course,  a  more  flexible  model  than  that  of  eq.  (12)  can  be  built  by 
introducing  a  distribution  function  for  the  thickness  of  the  dielectric  coating,  including  the  possibility  of  zero  thickness.  Such 
a  model  would  allow  for  contact  among  (only  a  fraction  of)  the  metal  spheres,  so  that  an  appropriate  choice  of  the  distribution 
function  should  ensure  a  smoother  transition  from  the  results  of  the  Bruggcman  model  and  those  of  the  Maxwell-Gamctt  model. 

The  results  we  discussed  above  show  that  our  Bruggeman-Mie  approach  yields  meaningful  results  even  when  applied  to 
particles  with  x  n  0.1  and  in  this  respect  proves  superior  to  the  customary  approaches,  eqs.  (6W8).  In  turn,  our  attempt  to 
introduce  some  correlation  among  the  metal  particles,  eqs.  (12),  proved  too  effective,  perhaps.  As  remarked  above,  the  amount 
of  correlation  w«  introduced  into  our  models  is  so  large  that  the  results  we  reported  can  give  only  semiquantitative  information, 
i.e.  they  can  indicate  the  trend  one  has  to  expect  in  more  realistic  situations,  that  could  be  effectively  simulated  along  the  lines 
outlined  above.  As  thev  stand,  the  models  we  presented  in  this  paper  can  give,  in  our  opinion,  useful  indications  on  the  dielectric 
behavior  of  the  composite  media. 
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Abstract 

Some  enhancements  to  the  Bruggeman  mixing  rule  designed  to  deal  with  com¬ 
posites  made  of  panicles  of  relatively  large  size,  possibly  correlated  in  their  distribu¬ 
tion.  arc  proposed.  The  resulting  implementations  arc  applied  to  a  composite  made 
of  metal  and  of  diclccric  spheres  cither  unconrclaicd  or  subject  to  a  correlation  of 
aggregation  among  the  metal  particles,  or  to  a  cojrTclaiion  of  exclusion  that  prevents 
contact  among  the  metal  spheres.  Our  calculations  show  that  our  enhanced  schemes 
work  well  up  to  sizes  ~  0 . 1  and  that  both  kinds  of  correlation  above  may  have  very 
strong  effects  on  the  transition  of  the  composite  from  the  dielectric  to  the  metallic 
behavior. 
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1  Introduction 


In  the  last  few  years  there  has  been  an  increase  of  interest  in  the  propagation  of  electromag¬ 
netic  waves  through  composite  media  in  view  of  the  practical  relavance  of  several  such 
materials  [  1 ,2,3].  A  possible  way  to  deal  with  this  propagation  problem  is  to  think  of  the 
macroscopic  properties  of  aj  composite  as  a  space  average  of  its  local  properties  and  conse¬ 
quently  to  substitute  the  actual  medium  with  an  effective  one.  In  practice  this  is  achieved, 
according  to  the  effective  medium  theory,  by  building  a  homogeneous  medium  capable  of 
fitting  the  average  propagation  properties  of  the  composite  in  such  a  way  that,  if  this  latter 
is  embedded  into  the  effective  medium,  no  scattering  is  experienced  {4].  There  is  no  uni¬ 
versal  recipe  to  build  the  dielectric  constant  of  the  effective  medium,  since  the  composites, 
by  their  very  nature,  may  liave  the  most  diverse  constitutions  and  topologies.  In  fact,  the 
properties  of  a  composite  depend  on  the  nature,  the  shape  and  the  size  as  well  as  on  the 
concentradon  and  on  the  possible  correlation  in  the  distribution  of  the  panicles.  A  partic¬ 
ularly  simple  model  composite  medium,  that  in  the  framework  of  the  effective  medium 
theory  leads  to  the  so  called  Bruggeman  mixing  rule,  is  composed  of  spherical  panicles 
of  very  small  size  parameter  with  totally  uncorrelated  space  distribution.  If  into  the  model 
above  is  forced  enough  correlation  as  to  ensure  that  the  panicles  are  well  separated  and 
dispersed  within  a  homogeneous  matrix  [2],  one  is  lead,  insted,  to  die  Maxwell-Gamett 
mixing  rule.  Of  course,  apart  from  such  an  extreme  treatment  of  the  correlation,  which 
leads  to  so  different  topologies,  the  approximations  of  the  two  models  above  are  rather 
severe.  In  actual  composite  media,  indeed,  the  size  parameter  of  the  particles  is  often  not 
very  small  and  the  panicleis  themselves  may  have  a  definitely  nonspherical  shape.  In  this 
respect  the  possibility  of  aggregation  of  several  spherical  panicles  should  also  be  consid¬ 
ered  as  this  phenomenon  produces  effectively  nonspherical  scatteren.  Even  by  assuming 
that  all  the  panicles  are  actually  spherical,  their  radii  are  seldom  identical  but  are  in  general 
distributed  within  a  more  dr  less  wide  range.  Of  course  some  work  has  been  done  to  over¬ 
come  at  least  some  of  these  approximations:  for  instance,  a  dynamical  extension  of  the 
Bruggeman  rule  applicable  to  size  parameters  up  to  0 .01  has  been  proposed  by  Stroud  and 
Pan  [4],  while  a  similar  extension  of  the  Maxwell-Gamett  rule  and  a  distribution  function 
for  the  radii  have  been  introduced  by  Chylek  and  Srivastava  [1]. 

In  this  paper  we  introduce  a  modified  Bruggeman  model  to  be  applied  when  the  size 
parameter  of  the  panicles  is  extended  up  to  0  1 .  Our  extension  is  effected  by  describing  the 
particles  through  the  exact  Mie  amplitudes  (5,6]  and  including  up  to  the  quadrupole  terms 
for  the  biggest  spheres.  We  tentatively  considered  also  a  medium  composed  of  panicles 
with  size  parameter  up  to  ~  l  .0  by  including  terms  up  to  the  2* -poles  to  ensure  the 
convergency  o'  our  calculations.  Nevertheless,  we  do  not  regard  the  results  we  obtained 
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for  this  latter  medium  as  reliable  although  they  come  out  as  a  natural  and  smooth  evolution 
of  the  results  for  spheres  with  smaller  size.  Indeed,  we  think  that  any  attempt  to  describe 
trough  the  effective  medium  theory  the  properties  of  a  medium  too  coarsely  grained  with 
respect  to  the  wavelength  of  the  propagating  radiation  is  bound  to  fail. 

We  apply  our  modified  model  to  study  also  the  effect  of  the  interparticle  correlation 
_on  the  behavior  of  the  percolation  of  a  composite  made  of  metal  and  of  dielectric  parti¬ 
cles.  The  first  kind  of  correlation  we  forced  onto  the  composite  medium  is  a  correlation 
of  exclusion  that  we  introduced  within  the  framework  of  our  Bruggeman-Mie  model  by 
putting  around,  each  metal  sphere  a  very  thin  dielectric  layer.  Then,  by  increasing  the 
thickness  of  the  layer  we  arc  able  to  follow  by  our  calculations  the  rather  sharp  decrease 
of  the  absorption  that  comes  with  the  implied  transformation  of  the  topology  from  that  of 
the  Bruggeman  model  to  that  of  the  Maxwell-Gamett  one. 

When  the  metal  particles  of  the  composite  are  allowed  to  aggregate  two  other  kinds  of 
correlation  set  on:  the  one  that  produces  the  aggregation  itself  and  the  possible  correlation 
in  the  orientations  of  the  resulting  anisotropic  particles.  We  assess  how  the  aggregation 
affects  the  dielectric  behavior  of  the  composite,  at  least  for  the  simple  case  of  binary  aggre¬ 
gation,  by  performing  some  calculations  through  an  approach  we  developed  in  the  last  few 
years.  In  fact,  that  approach  describes  the  scattering  by  aggregates  of  spherical  particles 
[71  through  multipole  amplitudes  analogous  to  those  of  Mie  as  well  as  the  macroscopic 
optical  constants  of  a  low-density  dispersion  of  such  aggregates  either  when  they  are  ori¬ 
ented  at  random  or  alike  [81.  Our  results  for  the  metal-dielectric  composite  referred  to 
above  show  that  the  aggregation  shifts  the  percolation  towards  lower  concentrations  of 
the  metallic  component  when  the  aggregates  are  randomly  oriented;  when  the  aggregates 
are  all  oriented  alike  larger  shifts  occur  with  their  magnitude  depending  on  the  direction 
of  incidence  and  on  the  polarization  of  the  propagating  wave. 

2  Theory 

The  effective  medium  theory  defines  etff  as  the  dielectric  constant  of  the  homogeneous 
medium  that  can  be  substituted  to  the  real  medium  without  affecting  the  propagation  of 
the  electromagnetic  waves.  This  amounts  to  say  that  ecff  should  satisfy  the  null  scattering 
condition ,  i.e.  it  is  the  solution  to  the  equation 

=  0  (1) 
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where  e,  is  the  polarization  vector  of  the  plane  wave  and  is  the  normalized  forward- 
scattering  amplitude  of  any  of  the  N<  identical  particles  of  the  t'-th  kind  when  immersed 
into  the  homogeneous  effective  medium  of  dielectric  constant  «efr. 

According  to  Mie,  the  normalized  forward-scattering  amplitude  of  a  homogeneous 
spherical  particle  with  refractive  index  n  *  e1/2  embedded  in  a  homogeneous  medium  of 
refractive  index  n„,  —  in  the  present  case  =  tv®  =  ('J(2  —  is  given  by 
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where  the  polarization  index  rj  has  been  dropped  for  the  sake  of  simplicity,  k  »  w/c  and 
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In  eqs.  (2)  and  (3)  the  ’s  are  the  muldpole  amplitudes  of  the  incident  plane  wave,  the 
’s  are  those  of  the  scattered  field  and 
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x  *  fcp  is  the  size  parameter  of  the  panicle  with  radius  p  and  p  is  a  parity  index  that  dis¬ 
tinguishes  the  magnetic  multipoles  ( p  -  1)  from  the  electric  ones  (p  =  2).  Accordingly, 
R^  refers  to  the  magnetic  and  Rf2)  to  the  electric  2*-pole.  When  |y«i|  <  1  eq.  (4)  can 
be  expanded  in  terms  of  ye„  and  one  gets 
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while  the  R(tp) 's  for  i  >  1  are  smaller.  As  a  consequence,  if  one  lets  the  size  parameter  x, 
of  the  particles  involved  in  eq.  (1)  tend  to  zero,  to  obtain  convergent  values  for  the  /,’ s, 
only  the  R\2)’s  as  given  in  eq.  (5)  need  be  considered;  then  eq.  (1)  turns  into  the  well 
known  Bruggeman  mixing  rule  [1] 
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where  the  x{  above  is  the  size  parameter  of  any  of  the  /*/,  spheres  of  volume  t \  and  dielectric 
constant  e,  and  VJ  =  Nxvt/V  is  the  fraction  of  the  total  volume  of  the  composite  that  such 
spheres  fill.  In  order  to  get  a  dynamical  extension  of  the  Bruggeman  mixing  rule,  Stroud 
and  Pan  [4]  do  not  neglect  the  R\l)  term  and  thus  include  into  the  expansion  of  eq.  (2) 
also  the  magnetic  dipole  term.  Since  they  take  R\l)  and  R\2)  in  the  form  of  eq.  (5),  their 
mixing  rule  reads 
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and  appears  as  an  extension  of  eq.  (6)  to  rather  larger,  although  still  small,  sizes.  The 
mixing  rule  of  Chylek  and  Srivastava  [1]  lies  on  the  very  same  approximations  but  these 
authors  include  the  size  dependent-terms  above  only  when  e  is  very  large  —  for  instance, 
this  is  the  case  of  metal  spheres.  They  also  include  a  size -distribution  function  X(p)  and 
write 
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Obviously,  the  size-dependent  terms  in  eq.  (8)  assume  the  same  form  they  have  in  eq.  (7) 
when  k;(p)  =  N,6(pj  —  p). 

Since  the  only  common  implementations  of  eq.  (1)  are  the  zero-size  approximation, 
eq.  (6),  and  the  small-size  approximation,  eqs.  (7)  and  (8),  we  determined  to  test  the 
reliability  of  the  null  scattering  condition,  eq.  (1),  when  applied  to  larger  particles.  In 
practice,  all  the  calculations  we  will  present  in  the  next  section  were  performed  by  carrying 
to  full  convergency  the  exact  Mie  expansions  to  the  effect  of  ensuring  the  inclusion  of  all 
significant  contributions  to  the  scattering  amplitude,  eq.  (2).  In  this  way  we  overcame  all 
the  limitations  due  to  the  small  size  of  the  panicles  but  those  that  prevent  the  possibility 
of  treating  the  whole  medium  through  a  dielectric  description. 

In  section  1  we  mentioned  the  special  kind  of  correlation  that  sets  on  into  the  space 
distribution  when  the  panicles  undergo  aggregation  phenomena.  Here  we  will  briefly 
describe  how  these  cases  can,  in  principle,  be  treated.  Indeed,  we  presented  elsewhere 
[7]  a  suitable  approach  to  the  scattering  properties  of  an  aggregation  of  spheres  (cluster) 
designed  to  yield  the  forward  scattering  amplitude  of  the  aggregate  as  a  whole  still  in  the 
form  of  eq.  (2)  but  for  that 


4l‘  =  -EEC'^’  (9) 
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where  the  quantities  defined  in  ref.  [9],  are  analogous  tothe  R\p)' s,  eq.  (3),  for 

the  single  spheres  and  take  account  of  the  specific  geometry  of  the  aggregation.  If,  as  it 
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often  occurs,  we  have  to  deal  with  a  dispersion  of  randomly  oriented  identical  aggregates 
and  /( 8)  is  the  scattering  amplitude  of  any  one  of  them  with  orientation  8,  we  can  sum 
over  all  the  orientations  and  get 

where  /  on  the  right  hand  side  is  still  in  the  form  of  eq.  (2)  but  in  this  case 

- -yfi t  E  [e  s£K.-<*>1  vg>.  < io) 

Thanks  to  the  structure  of  the  summation  in  eq.  (10),  the  resulting  ^4’s  are  independent  of 
the  particular  orientation,  8o ,  we  chose  to  perform  the  calculation.  The  quantities  S  in  eqs. 
(9)  and  (10)  are  the  elements  of  a  matrix  that  must  be  obtained  by  inversion  of  a  square 
matrix  of  order  2  X.L{  L  +  2) ,  defined  in  ref.  [8],  where  *1  is  the  number  of  the  spheres  in 
the  aggregate  and  L  is  the  value  of  £  at  which  the  sums  are  truncated.  Considerations  on 
the  symmetry  of  the  aggregate,  if  any,  can  be  used  to  factorize  the  matrix  to  be  inverted 
[10,11];  in  this  respect,  the  most  convenient  case  occurs  when  the  centers  of  the  spheres 
lie  on  a  straight  line,  for  the  maximum  order  of  the  matrices  to  be  inverted  is  2^L  only. 
In  the  next  section  we  will  deal  with  binary  aggregates  with  a  truncation  £  appropriate  to 
the  size  of  the  spheres  so  that  L  <  4  has  been  sufficient  in  all  cases  and  the  computational 
effort  has  been  therefore  rather  low. 

As  mentioned  in  the  preceding  section,  we  also  considered  a  model  composite  medium 
which,  in  a  sense,  is  a  hybrid  between  that  of  Bruggeman  and  that  of  Maxwcll-Gamett; 
namely  we  calculated  ce n  for  a  composite  made  of  homogeneous  dielectric  spheres  with 
«  =  and  of  spheres  with  a  metal  core,  of  dielectric  constant  e u .  coated  by  a  thin  dielec¬ 
tric  layer.  The  scattering  properties  of  the  coated  spheres  were  not  described  through  an 
approach  based  on  the  boundary  conditions  at  the  separation  surface  but  rather  by  inter¬ 
posing  between  the  core  and  the  coating  a  very  thin  transition  layer  [9]  through  which 

e  =  e(  r)  =*  tu  +  ( 3  a2  -  2  s3 )  A 


with 


A  =  «o  —  and 


r  -  Pm 
Pd  -  Pm' 


where  py  is  the  radius  of  the  metal  core  and  pn  the  inner  radius  of  the  dielectric  coating.  In 
this  way  we  achieved  a  smooth  transition  between  the  two  materials  as  the  refractive  index 
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n  =  n(  r)  is  a  regular  function  of  r  everywhere  within  the  coated  spheres.  Consequently 
eq.  (4)  must  be  substituted  by  [12,6] 

(t)  C'/’u,  -(14. 

R‘  1 

“where  £^p>  s*  ( x)  are  now  calculated  by  numerical  integration  of  the  radial  equations 
d’G*'1  2  dnd{3$:) ,  f,  ««+ Dl  „ 

~W~  -  n  dJ-dT8’2  +  l"  “  ~~p~J  ‘ 

with  £  =  kr.  Of  course,  when  the  spheres  are  homogeneous,  «o  3  andeq.  (11)  reduces 
to  eq.  (4).  As  a  matter  of  fact,  all  the  calculations  we  will  present  in  the  next  section  were 
performed  by  using  eq.  (11),  or  eq.  (4)  when  appropriate,  with  no  approximation.  It 
is  worth  noticing  that  our  previous  calculations  [9]  proved  that  the  introduction  of  the 
transition  layer,  provided  it  is  thin  enough,  has  no  influence  on  the  results. 


3  Results  and  discussion 


All  the  composite  media  we  dealt  with  in  our  calculations  are  two-component  ones;  the 
first  component  is  a  dielectric  with  eD  =  1  and  the  second  is  a  metal  with 


3  1 


1 

i/(  u  +  »7)  ’ 


(12) 


where  v  =  u/uip,  7  =  1  /(  tup),  viz.  the  dielectric  material  is  the  vacuum  while  the  metal¬ 
lic  component  is  described  by  the  free-electron  Drude  constant  [13].  This  choice  allows 
us  an  easy  comparison  with  the  preceding  results  of  Stroud  and  Pan  [4]  and  of  Chylek  and 
Srivastava  [1].  Anyway,  we  think  that  a  greater  sophistication  of  the  description  of  the 
metallic  component  is  of  no  usefulness  within  a  framework  substantially  as  naive  as  the 
effective  medium  one. 

If  we  assume  the  particles  to  be  spherical,  eq.  (1)  becomes 


Ktr/w(«eff)  +  (Pw/Pd)3(  1  -  -  0,  (13) 


where  p y  is  the  radius  of  the  metal  spheres  and  po  the  radius  of  the  dielectric  ones,  and 
the  f's  are  given  by  eqs.  (2),  (3)  and  (4). 
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When  the  metal  spheres  are  coated  by  a  dielectric  layer  we  have 

Vufuc(ttti)  +  (Pw/Po)3!  1  —  Vu(  1  +  Ci)]/o(«cff)  =0,  (14) 

where  fuc  refers  to  a  sphere  with  a  metal  core  of  radius  pu  and  a  dielectric  coating  layer 
with  external  radius 

Pi  *  (1  +  cL)ippu, 

and  is  thus  given  by  eqs.  (2),  (3)  and  (11).  The  quantity  ci  is  defined  as 

cl  =  vl/vu 

where  vl  is  the  volume  of  the  dielectric  layer:  in  other  words  a  gives  the  volume  of  the 
dielectric  coating  around  each  metal  sphere  as  a  fraction  of  the  volume,  vu,  of  the  metal 
spheres  themselves.  Of  course 

5s  =  .  ”t  0  ^  C[,  <[  Cimui  ~  '  1, 

1  +  cL  Vu 

for  a  given  c*,  and  for  a  given  Vu,  respectively,  so  that  when  ct  =*  0,  i.e.  when  the 
spheres  have  no  coating,  we  get  again  the  topology  of  the  Bruggeman  model,  whereas, 
when  Vu  -  Vum»t  or  ct  =  ctm,»  we  get  the  Maxwell-Gamett  topology. 

When  the  spheres  aggregate  in  pairs  we  have 

Vufcit'")  +  2 (pu/PoVl  1  -  Vwl/D(fcff)  -  0,  (15) 

where  fc  is  the  contribution  of  one  aggregate  averaged  over  the  orientations  as  described 
in  the  preceding  section. 

It  is  perhaps  worth  noticing  that  eqs.  (13MI5),  unlike  the  two-component  implemen¬ 
tations  of  eqs.  (6)-(8),  are  transcendent  equations  in  eta,  and  thus  require  some  further 
caution  in  seeking  for  their  solutions. 

We  calculated  c€ff  according  to  eqs.  (13),  (14)  or  (15)  as  a  function  of  Vu  for  some 
values  of  xu  =  kpu  and  of  ir,  7  =  0  01  in  all  cases.  In  particular  we  report  in  figs.  1 
and  2  Imfetffl  as  calculated  from  eq..  (13),  and  from  the  two-component  implementations 
of  eq.  (6),  of  eq.  (7)  with  po  -  Pu,  and  of  eq.  (8)  with  a  <5-Iikc  distribution  of  the  radii. 
Our  main  purpose  is  to  show  the  most  important  features  of  the  absorption  of  a  composite 
medium  at  frequencies  for  which  the  metal  component  is  by  itself  strongly  absorbing. 
Therefore,  we  chose  to  plot  our  results  for  v  =  0  01 ,  Fig  1,  and  fori/  =  0.1,  Fig.  2,  since, 
for  many  metals,  these  frequencies  occur  in  the  far  and  in  the  near  infrared,  respectively. 
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It  is  worth  noticing  that  both  our  eqs.  (13),  (14)  and  (15)  and  their  counterpart  according 
to  Stroud  and  Pan,  eq.  (7),  show  a  dependence  on  the  size  of  the  dielectric  spheres.  When 
one  assumes  tu  =  1 ,  i.e.  when  the  dielectric  is  the  vacuum,  the  problem  arises  whether  its 
more  or  less  minute  subdivision  into  spherical  cells  may  affect  the  results.  As  a  matter  of 
fact,  we  took  xo  =  kpo  =  10  "4  but  our  calculations  confirm  that  the  size  of  the  vacuum 
.  spheres  has  very  small  influence  either  when  their  size  is  reduced  to  zero  according  to 
"  Chylek  and  Srivastava  or  when  it  is  set  equal  to  the  size  of  the  metal  spheres  according  to 
Stroud  and  Pan.  As  expected,  our  results  are  nearly  independent  of  the  size  of  the  vacuum 
spheres,  indeed,  for  all  the  sizes  of  the  metal  spheres  at  which  the  implementations  of  the 
effective  medium  theory  we  deal  with  can  be  assumed  to  work  correctly.  Of  course  this 
is  true  as  far  as  ctg  does  not  become  so  large  within  the  range  of  Vy  as  to  enhance  the 
size-dependent  contributions  from  the  dielectric  speres.  In  practice,  the  range  of  Vy  both 
in  Fig.  1  and  in  Fig.  2  is  not  too  large  and  yet  it  extends  enough  to  show  the  influence 
of  the  size  and  of  the  frequency  on  the  transition  from  a  dielectric  behavior  to  a  metallic 
one.  In  our  opinion  xy  -  0.1  is  an  appropriate  upper  limit  for  the  applicability  of  eq. 
(13),  viz.  it  is  a  practical  limit  to  the  size  of  the  grains  of  a  composite  medium  in  order 
that  this  latter,  from  the  point  of  view  of  the  dielectric  properties,  could  be  substituted  by 
a  homogeneous  medium.  Anyway,  on  the  basis  of  our  calculations  we  conclude  that,  for 
particles  of  size  greater  that  0.1  but  rather  smaller  than  l,theformofIm[eeff]  ffomeq.  (13) 
is  still  meaningful  in  the  small-  Vy  region.  In  this  connection  we  recall  that  eq.  (6)  should 
never  be  applied  when  x  >  0 .01  and  that  eqs.  (7)  and  (8)  cannot  be  thought  of  as  accurate 
when  xy  ~  0  1 ,  because  these  equations  come  out  from  a  series  expansion,  eq.(5),  of  the 
exact  Mie  quantities,  eq.  (4).  We  notice  that  in  Fig.  2  the  curves  for  xy  =  0.1  from  eqs. 
(7)  and  (8)  coincide  with  each  other  and  are  substantially  different  from  that  coming  from 
eq.  (13),  while  they  do  not  appear  at  all  in  Fig.  1  for  the  approximations  on  which  eqs. 
(7)  and  (8)  are  b->sed  result  in  solutions  that  are  physically  unacceptable.  In  fact,  these 
equations,  that  an  algebraic  equations  of  the  third  degree,  yield  values  of  Im[  «rff]  that  are 
either  negative  or  everywhere  vanishingly  small  or  very  large  even  in  the  low  Vy  region. 

In  Figs.  3and4wepIotIm(<eff]  fort/  » ‘V01  and  for  several  values  of  a  as  calculated 
according  to  eq.  (14);  analogous  results  were  obtained  for  u  ■  0.1.  It  is  quite  evident  that 
the  correlation  of  exclusion  produced  by  the  dielectric  coating  around  the  metal  spheres 
has  very  strong  effects  even  at  very  low  values  of  c£.  Actually,  the  model  described  in 
eq.  (14)  prevents  any  contact  among  the  metal  particles  and  is  therefore  topologically 
akin  to  the  Maxwell-Gamett  model.  As  a  consequence,  it  is  not  surprising  that,  when  ci 
increases,  our  results  tend  very  fast  to  those  that  are  typical  of  the  latter  model.  Figs.  3 
and  4  show  the  results  for  id  =  10  .  At  this  value  of  xg  the  contribution  of  the  vacuum 

spheres  to  the  scattering  amplitude  is  very  well  approximated  by  the  first  of  of  eqs.  (5) 
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and  such  a  small  size  suggests  the  idea  of  a  dielectric  medium  which  pervades  and  fills  all 
the  space  around  the  coated  spheres.  Even  in  this  case  any  size  of  the  dielectric  spheres 
up  to  xd  ■  xu  proved  to  be  incapable  of  affecting  the  results  of  our  calculations  from  eq. 
(14).  Of  course,  a  more  flexible  model  than  that  of  eq.  (14)  can  be  built  by  introducing 
a  distribution  function  for  the  thickness  of  the  dielectric  coating,  including  the  possibility 
of  zero  thickness.  Such  a  model  would  prevent  the  contact  among  (only  a  fraction  of)  the 
metal  spheres,  so  that  an  appropriate  choice  of  the  distribution  function  should  ensure  a 
smoother  transition  from  the  results  of  the  Bruggemah  model  and  those  of  the  Maxwell- 
Garnett  model. 

Figures  5  and  6,  where  Im[eeff]  is  plotted  from  eq.  (15),  for  u  *  0.01,  show  the 
effect  of  a  complete  binary  aggregation  of  the  metal  particles  on  the  transition  of  a  com¬ 
posite  medium  to  the  metallic  behavior,  analogous  results  were  obtained  for  v  *0  1. 
This  kind  of  effect  proves  a  quite  significant  one  but  it  is  very  strongly  increased  when 
to  the  correlation  of  aggregation  one  also  adds  the  request  that  all  the  binary  aggregates 
be'  likely  oriented.  Even  in  this  case  the  value  of  zD  proves  to  be  irrelevant,  although  the 
conceptually  most  satisfactory  choice  is  that  of  very  small  size  (in  practice  iD  =  10  -4). 
Both  Figures  5  and  6  show  that  the  percolation  undergoes  a  noticeable  shift  towards  higher 
values  of  Vy  when  the  common  orientation  of  the  aggregates  is  chosen  so  as  their  axes 
are  parallel  to  the  wavevcctor  of  the  incident  wave.  This  phenomenon  is  easily  under¬ 
standable  if  one  thinks  of  the  arrangement  above  as  the  one  that  provides  the  maximum 
of  metal-free  paths  to  the  propagating  radiation.  On  the  other  hand,  the  maximum  in  the 
extinction  power  of  the  medium  that  occurs  when  the  the  axes  of  the  aggregates  are  all 
parallel  to  the  electric  field  can  also  be  attributed  to  the  arrangement,  as  this  choice  of  the 
orientation  maximizes  the  absorption  of  the  individual  aggregates.  Of  course,  less  strong 
effects  should  be  expected  if,  more  realistically,  only  a  reasonable  fraction  of  the  metal 
particles  were  allowed  to  aggregate. 

The  results  we  discussed  above  show  that  our  Bruggeman-Mie  approach  yields  mean¬ 
ingful  results  even  when  applied  to  pdrticles  with  x  ~  0 . 1  and  in  this  respect  proves 
superior  to  the  customary  approaches,  eqs.  (6M8).  In  turn,  our  attempt  to  introduce  some 
correlation  among  the  metal  particles,  eqs.  (14)  and  (15),  proved  too  effective,  perhaps. 
As  remarked  above,  the  amount  of  correlation  we  introduced  into  our  models  is  so  large 
that  the  results  we  reported  can  give  only  semiquantitative  information,  i.e.  they  can  in¬ 
dicate  the  trend  one  has  to  expect  in  the  more  realistic  situations  that  could  be  effectively 
simulated  along  the  lines  outlined  above.  As  they  stand,  the  models  we  presented  in  this 
paper  can  give,  in  our  opinion,  useful  indications  on  the  dielectric  behavior  of  the  com¬ 
posite  media. 
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Figure  captions 

Figure  1  :  Im[ceff]  vs.  Vw,  from  eq.  (13)  for  v  =  0.01  andz d  3  10~4.  The  results  from 
eqs.  (6>— (8)  are  also  reported  for  the  sake  of  comparison.  The  label  B  denotes  the 
curve  from  eq.  (6);  as,  on  this  scale,  several  curves  are  undistinguishable  from  each 
other  the  labels  give  xu  and  the  number  of  the  equations  that  generate  the  curve. 

Figure  2  :  Same  as  Fig.  1  but  for  v  -  0 . 1 . 

Figure  3  :  Im[c<ff]  vs.  Vy,  from  eq.  (14),  calculated  for  the  values  of  ci  that  label  the 
curves.  =  0 .01,  xp  3  10-4  and  i>  =  0 .01. 

Figure  4  :  Same  as  Fig.  3  but  for  x  u  3  0  1 . 

Figures  :  Im[eefr]  vs.  Vw,  from  eq.  (15),  for  binary  aggregates  with  xu  »'  0  01, 
x  a  =  10  “4  and  v  =  0 .01 .  The  uniabelled  solid  curve  refers  to  randomly  oriented 
aggregates;  that  labelled  ||  refers  to  aggregates  all  oriented  with  their  axes  parallel 
to  the  wave  vector  and  the  one  labelled  ±||  refers  to  aggregates  with  their  axes  par¬ 
allel  to  the  electric  field  of  the  wave  which  is  assumed  to  be  linearly  polarized.  The 
dashed  curve  repons,  for  the  sake  of  comparison,  the  results  from  eq.  (13). 

Figure  6  :  Same  as  Fig.  5  but  for  x.w  =0.1. 
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Abst net :  The  most  significant  results  we  obtained  on  the  propagation  of 
electromagnetic  waves  through  non  homogeneous  media  are  summarized. 


1 .  Introduction 

A  few  years  age  we  devised  a  model  scatterer  suitable  for  a  systematic 
study  of  the  optical  properties  of  aerosols  and  in  particular  of  the  effects 
of  the  anisotropy  and  of  the  possible  aggregation  of  the  constituent  parti¬ 
cles.  Indeed,  our  model  scatterer  is  built  as  a  cluster  of  spherical  objects 
of  arbitrary  radii  and  retractive  indexes;  the  spheres  need  not  be  homogene¬ 
ous  md  their  relative  positions  are  arbitrary  to  a  large  extent.  The  scat¬ 
tering  st  electromagnetic  waves  by  such  composite  scatterers  is  dealt  with 
through  the  use  of  multipolo  expansions  to  describe  both  the  incident  and 
the  scattered  field  as  well  as  the  field  within  the  component  spheres  (l,  2). 
The  resulting  forward  scattering  amplitude  has  very  simple  t rans f o rmat ion 
properties  under  rotation  of  the  coordinate  axes.  7hi3  allows  us  to  calculate 
very  easily  the  macroscopic  optical  properties  of  a  low-density  dispersion 
:f  clusters  even  when  they  are  randomly  oriented.  In  fact,  we  are  able  e.g. 
to  3tudy  the  changes  of  the  spectrum  of  a  dispersion  of  spherical  particles, 
both  homogeneous  and  r.on-hcmoqeneou3  13,4],  when  they  undergo  various  stages 
of  aggregation  and  even  when  they  change  their  mode  of  aggregation  in,  analo¬ 
gy  to  what  happens  ir.  chemical  reactions  (5). 

In  the  iu3t  year  we  started  an  attempt  to  extend  our  study  to  intermedi¬ 
ate  and  high-den3ity  dispers ions  through  a  modified  version  of  the  effective 
medium  theories  i  6  j  .  Essentially  the  effective  medium  theories  substitute  the 
actual  dispersion  with  an  effective  medium  whose  dielectric  properties  are 
calculated  so  as  to  compensate  exactly  for  the  scattering  produced  by  the 
particles.  Cur  particular  work  in  this  field  consists  in  treating  a  disper¬ 
sion  of  metal  spheres  with  size-parameter  bigger  than  that  usually  dealt 
with  within  the  Bruggemann  scheme  [7]  and  yet  sufficiently  small  to  ensure 
the  validity  of  the  dielectric  description  of  the  medium.  Accordingly  we  used 
the  full  Mie  expansion  and  as  a  further  improvement  included  some  kind  c£ 
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correlation  «.g.  that  produced  by  aggregation  phenomena  or  by  an  exclusion 
hole  around  each  sphere.  This  amounts  to  using  a  mixed  scheme  midway  between 
that  of  Sruggemann  and  that  of  Maxwell-Garnett .  Some  of  our  preliminary  re¬ 
sults  will  be  illustrated  in  Section  3. 


2 .  Theory 

All  the  information'  on  the  propagation  of  electromagnetic  waves  through 
a  low-density  dispersion  of  scatterers  is  contained  in  the  matrix  of  the  re¬ 
fractive  index  18] 


A/-, 


In 

VV 


11) 


where  the  sum  runs  over  all  the  scatterers  within  the  volume  V  and 

/..w, u * 121 


In  eq.  (2)  is  the  polarization  vector  and  fVT^  the  forward-scattering  am¬ 
plitude  of  the  v-th  scatterer;  in  turn  the  *s  are  the  multipole  am¬ 
plitudes  of  the  incident  plane-wave  field  and  che  a<p,LM's  are  their  counter¬ 
parts  for  the  scattered  field.  The  superscript  p  is  a  parity  index  which  dis¬ 
tinguishes  the  magnetic  muitipole  fields  (p-1)  from  the  electric  ones  (p*2) 
(9| .  When  the  scatterer  is  a  cluster  of  spheres  the  muitipole  amplitudes  of 
the  scattered  field  are 


j>  _  v-*  r/O'V) 

- 2L*  U.LU^L’U' 

yLU- 


where  the  S' s  are  defined  in  ref. (3).  The  most  important  feature  of  tho  S's 
is  that  under  rotation  of  the  coordinate  axes  they  transform  according  to  the 
representations  of  the  full  rotation  group.  As  a  consequence,  once  the  S's 
are  known  for  a  cluster  of  given  orientation  they  are  immediately  known  for 
any  other  cluster  of  different  orientation,  and  the  sum  in  eq.(l)  can  be  per¬ 
formed  analytically. 

For  intermediate  and  high-d°nsity  dispersions,  the  procedure  outlined 
above  does  not  yield  reliable  results  because  the  multiple-scattering  pro¬ 
cesses  become  more  and  more  important.  A  possible  way  to  overcome  this  diffi¬ 
culty  is  to  resort  to  the  effective  medium  theories  such  as  the  well-known 
mixing  rules  of  Bruggeman,  and  of  Maxwell-Garnett'.  The  general  principle  on 
which  these  rules  are  based  was  outlined  in  sect.  1  and,  as  our  work  in  the 
field  is  still  in  progress,  we  defer  our  further  comments  to  the  next  section 
where  some  our  preliminary  results  will  be  presented. 
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Tig.  1  we  report  the  real  (la)  and  the  imaginary  part  (lb)  of  Njj  for  a 
low-density  dispersion  of  binary  clusters  both  oriented  alike  (solid  lines) 
and  at  random  (broken  lines)  versus  the  size  parameter,  x-kp,  of  the  constitu¬ 
ent  spheres,  which  are  homogeneous  with  n*1.3.  All  the  curves  are  normalized 
to  the  corresponding  quantities  of  a  dispersion  of  independent  spheres  with 
the  same  density  and  refractive  index  and  show  not  only  a  noticeable  depen¬ 
dence  on  the  orientation  but  also  that  for  a  cluster  of  given  geometry  there 
exist  a  value  of  the  size  at  which  this  dependence  disappear'.  Our  experience 
suggests  that  this  pseudo-spherical  behavior  is  bound  to  appear  for  some  value 
of  the  size  when  the  clusters  have  a  definite  symmetry. 


As  anticipated  in  3ect  .  1,  we  also  dealt  with  dispersions  of  layered 
spheres  which  we  treated  as  explained  ; n  ref. [3).  In  short,  w«  kept  constant 
the  dielectric  function  within  each  layer  and  between  each  pair  of  them  inter¬ 
posed  a  transition  layer,  as  th\n  as  possible,  within  which  the  radial  depen¬ 
dence  of  the  dielectric  function  is  ta*en  so  as  to  ensure  the  continuity  of  C 
and  of  its  radial  derivative,  w*?  considered  in  this  way  spheres  with  a  diffuse 
surface,  metal  sphere?  with  a  dielectric  or  metallic  coating  and  dielectric 
spheres  wi’th  a  metallic  coating.  The  dielectric  function  of  the  metal  was  de.- 
scribed  by  a  f ree-electron  Drude  function,  while  that  of  the  dielectric  was 
described  by  a  damped  oscillator  function.  The  most  striking  result  is  shown 
in  fig.  2  which  refers  to  spheres  of  M.q 0  coated  with  Al  with  a  total  radius  of 
50A.  It  is  quite  evident  that  the  absorption  peaks  of  the  dispersion  of  single 
spheres  (solid  lines)  shift  according  to  the  thickness  of  the  coating:  thus 
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this  feature  can  be,  in  principle, 
used  to  obtain  a  selective  absorp¬ 
tion  at  any  frequency  within  a  given 
range.  Nevertheless  this  effect  can 
be  greatly  weakened  if  a  considera¬ 
ble  percentage  of  the  spheres  aggre¬ 
gate  (broken  lines)  . 

At  last  we  give  here  an  antici¬ 
pation  of  our  present  attempt  to  use 
at  best  the  potentialities  of  the 
effective  dielectric  constant  meth¬ 
ods  beyond  the  mixing  rule*  of 
Bruggeman  and  of  Maxwel 1-Garnett .  In 
fig.  3  we  present  our  results  for 
one  of  the  conceivable  intermediate 
models  between  that  of  Bruggeman  and 
that  of  Maxwell-Ga rnett .  This  model 
consists  of  an  admixture  of  metal 
and  dielectric  spheres  (these  latter 
with  dielectric  constant  e- 1  )  . 
Around  each  metal  sphere  we  put  a 
exclusion  hole  in<  the  form  of  a 
layer  of  variable  thickness  with 


Fig. 
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dielectric  constant  t-l.The  curves  in  fig.  3  are  labeled  by  the  percentage  of 
volume  occupied  by  the  exclusion  holes.  The  curve  labelled  B  refers  to  the 
conventional  Bruggeman  rule.  All  the  spheres  were  dealt  with  through  the  full 
Mie  theory.  The  main  result  of  our  Calculation  is  a  dramatic  shift  of  the 
per  eolation  threshold  as  a  consequence  of  the  correlation  of  exclusion  we 
forced  on  the  system. 
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ABSTRACT 

A  ee lection  of  clusters  compoeed  of  20  spheres,  simnlnting  irregularly  shaped 
particles  often  observed  in  powders,  are  used  to  model  a  low  density  monodispersed 
aerosol  for  which  the  absorption  coefficient  is  calculated.  Alto  calculations  will  be 
shown  for  small  clutters  (  compoeed  of  2  spheres)  havinf  the  tame  optica] 
characteristics  as  thoea  of  existing  experimental  data  with  which  comparison  evil]  be 
made. 


1.  Introduction. 

The  fact  that  aeroaol  particlae,  bat  for 
few  exception!,  are  not  spherical,  nor  can  tie 
described  as  spherical  in  consideration  of  their 
optical  properties,  is  well  known.  It  is  not 
surprising  that  tbs  optical  behaviour  of  aeroaol 
particles  as  calculated  using  Mis  Theory  is  eel* 
dom  consistent  with  the  experimentally 
observed  data,  sine*  ths  Mis  theory  is  based 
upon  the  sphericity  of  the  ecalterer.  It  la  how¬ 
ever  important  to  achieve  a  reasonable  confi¬ 
dence  in  solving  the  problem  of  scattering  by 
an  irregularly  shaped  particle.  The  problem  is 
not  easily  solved  and  special  difficulties  arise 
when  the  ecalterer  is  equal  to  or  larger  than 
the  wavelength  of  the  incident  light  and  when 
the  index  of  refraction  of  the  surrounding 
medium  differs  mors  than  for  a  small  amount. 
We  have  developed  a  formalism  which  allows 
ns  to  describe  an  irregularly  shaped  parti- 
,  and  the  calculation  of  the  properties 


of  snob  a  complex  sc  alter  tr  can  be  performed 
following  the  approach  outlined  in  previous 
works*"1!  Having  obtained  in  thin  manner  the 
optica]  properties  of  a  single  scsttrrer  then  the 
next  step  it  to  calculate  the  macroscopic  opts- 
cal  roast  an  la  for  the  propagation  through  n 
model  aerosol.  We  will  show  that  our  compu¬ 
tation  for  two  sphsres  agree  with  experimental 
results^,  than  w«  will  discuss  the  case  of  some 
complex  atructurei  which  are  types  found 
commonly  in  powdsrsl1*"1^ 

2.  Theory. 

We  refer  yon  to  other  papers  for  the 
details  of  the  computations  and  their  justifies- 
tiona,  ,  here  simply  accept  that  wt 

can  writ#  the  extinction  croea  section  of  a 
whole  cluster,  where  a  duster  is  the  geometri¬ 
cal  description  of  a  particle,  in  the  form 

<r,„  -  <  Vt|S|V  >  (1) 
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when  tk«  matrix  S  diptidi  oatjr  ox  tha 
XneMn  of  tka  dealer. 

Sa  1 1 1  x  )T|  xe***  (2) 


Tka  matrix  §  iaelmdae  all  tka  iaformuioa  oa 
tka  ac  at  taring  powar  of  tka  daatar  u  a  whole  ( 
ila  alamaata  depend  only  oa  tka  featuring 
powar  of  tka  apkaraa  of  tka  daatar  aad  oa  tka 
ralativa  poaitioaa  of  tkair  can  tan).  Ia  aq.(2), 
tka  matticaa  1  aad  T  wart  defined  prrrv 
ooatjr^.  Matrix  I  kaa  tka  rola  of  referring  tka 
amplitadaa  af  tka  maltipola  Saida  k  at  tar  ad  by 
aack  apkara  to  tka  orifia  of  tka  coordiaataa. 
Tka  matrix  T  of  tka  (bra: 


ju.il  ju.il 
ju.x  jux 


(S) 


wkick  mart  ba  cakalatad  aa  tka  tavaraa  af  tka 
matrix : 


M  - 


piCIJ-'+H  X 

I  '  El'T'+H 


(«). 


aad  wkick  ia  dalarmiaad  if  tAa  atraaara  a/  Ua 
ilula'  Tka ,  matrix  V  rapraaaata  tka 
amplitadaa  el  tka  maltipola  Saida  aad  ita  ala¬ 
maata  make  eq.(l)  dapandaat  oa  tka  diractioa 
of  iacidaaca.  Tka  matrix  H  daacribaa  ia  tka 
vicinity  el  tka  aorfaca  el  aack  apkara  tka  mak 
tipola  Saida  acattarad  by  tka  otkar  apkaraa  of 
tka  daatar.  Tka  matrix  X  rapraaaata  tka 
croon  tarma  of  tka  maltipola  Saida.  Aad  finally 
tka  matricaa  X,1|andX111  ara  tka  neutering 
powar  of  tka  aiagla  apkaraa  compoaiag  tka  daa¬ 
tar  (  idaatical  la  phyaical  moaning  to  tka 
tarma  A  aad  B  ia  Van  Da  Halat|L>*  ). 

I.  Xaralta  aad  Diacnaaloa 

Two  apkaraa  —  touching.  Tka  Snt  tat 
of  data  aril]  ba  that  of  tka  two  apkaraa,  which 
will  ba  comparad  with  tha  experimental  data 
of  Shuarmaa  aad  Wang'*'.  Tka  apharaa  era 
tooching,  tha  qaantitwa  which  wa  kava  tom- 
patad  ara  tha  complax  forward  acatlaring 
amplitada  |  S(0)  |  aa  a  function  of  tha  aiaa 
parameter  ha  *  2 n/\  el  tha  daatar.  Tha 
gaomatry  aad  all  tha  paramatara  wa  hava  aaad 
aad  that  wa  report  ara  tha  tuna  aa  tha  expnri- 


Bteata  for  wkick  wa  compara.  Following  we 
raprodnea  tka  axparimaatal  tr*»p: 


tka  following  grapk  ikowa  oar  ealcalntad  data 
(•olid  uat)  aad  tka  axparimaatal  dun  (dottad 
ttna). 


ka 


N  Spharaa.  Wa  will  tow  couaidar  dun-  . 
tan  com^ooad  of  20  apkaraa  in  ordtr  to  timn- 
iau  coma  irregularly  ahapad  aggregation!  oftan 
found  ia  powdery  In  the  Sgaroo  we  report 
tka  qaantitier  Twy/q  „  and  T „m1mh  w  f* 
aavaral  direction  of  incidence  u  a  fuaction  of  x 
tka  aiaa  parameter  ;  where  7„  ia  tha  abaorp- 
tioa  coefficient  of  a  diapenioa  of  apkaraa  of 
aiaa  parametr  y1>w201^  x<  whoa*  number  den¬ 
sity  tqaali  thu  of  the  diaperaion  of  duaten. 
Tha  quantity  7  it  tha  abtorptioo  coefficient  'or 
tha  daatar  randomly  oriented  and  7.  depend* 
oa  the  orientation  of  tha  cloatera  with  reaped 
to  tha  incident  plana— wave.  The  material  of 
tha  apharaa  compoaing  tha  daatar  in  all  cate* 
waa  chaaaa  to  ba  homoganaoaa  iaotropic, 
dielectric,  aad  nondiaparriva  arith  refractive 
index  wel.J  aad  finally  all  tha  apharaa  ia  tha 
doattn  are  identical. 
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The  clutter  in  defined  an  a  Strip  and  ita 
gtomatric  (aaturea  are  ahown  above,  ,1a  the 
graph  below  the  clutter  T  (broken  carve)  and 
r.  (tolid  linen)  veraua  x  lor  never al  directiona 
of  incidanca  for  a  dupe  re  ion  of  atxipa.  The 
incident  plane  wave  in  circularly  polariaed. 
The  angina  of  incidence  (teta,d)  ,  in  degreea, 
Labeling  the  corvee  an  related  to  the  com¬ 
mon  orientation  of  the  cluatere 


The  c hitter  it  defined  an  a  Empty  rectangle 
and  ita  geometric  feat  or  an  are  ahown  above.  In 
the  graph  below  the  dorter  V  (broken  curve) 
and  r.  (tolid  liaaa)  venae  *  far  aeveral  direc¬ 
tion  of  incidence  for  a  diaperoioa  of  empty 
ractaaglea.  The  incident  plane  wave  in  circa- 
larly  polariaed.  The  angler  of  incidence 
(leta,d)  ,  in  degreea,  labeling  the  Vm  curvet 
are  related  to  the  common  orientation  of  the 
duaUra. 


w 


f.li 


The  duiv 
feometric 
fraph  below 
T.  (solid  Urn 
of  iacideace  I 
incident  plai 
Th»  an | In  t 
labeling  th« 


mon  onrntatioi 


figure  I 

defined  a*  a  Rte Unfit  tad  ita 
area  aro  shows  above.  la  tbe 
[the  clutter  r  (brakea  carve)  aad 
a)  venae  X  for  several  directions 
ir  a  diapenioa  of  rectangles.  Tbe 
e  wave  ia  circslarly  polarised. 
F  iacideace  (lets,*)  ,  ia  degrees, 
'm  curves  are  related  to  the  com- 
«  of  the  clusters. 


Ifni 

The  duster  ia  defused  aa  a  Mriwerh  aad  tU 
feometric  feat  area  an  above  above.  Ia  the 
raph  below  the  chaster  T  (broker  carve)  aad 
r.  (a°W  liaea)  verves  X  fa r  several  directioaa 
of  iacideace  lot  a  diapenioa  of  art  weeks.  The 
iacidrat  plaae  ware  ia  ocularly  polarised. 
The  aailea  of  iacideace  (tata,b)  ,  ia  degress, 
labeliaf  the  P„  corvee  aae  related  to  the  com- 
moa  cereal  at  toe  of  the  dealers. 
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flgun  T  flgun  • 

Tkt  clatter  is  defined  a*  a  Raafinp  hie  and  it*  Tka  clutter  it  dtfiaed  M  a  Mngular  itrip  and 

geometric  baton*  an  thow*  abort.  Ia  iki  «*  geometric  baton*  art  thow.  abort,  la  tbt 

graph  below  tbt  cloatar  T  (broken  curve)  tad  P»pk  below  tbt  tloatar  T  (broken  cam)  and 

r-  (*obd  lint*)  nrtot  x  for  aartral  direction*  I*,  (to lid  Unta)  venua  x  for  tavaral  direction* 

of  incidence  for  a  diapenioa  of  roofing  tilt*.  °f  incidence  for  a  dijparaioa  of  angular  itripa. 

Th*  incidtal  plan*  war*  ia  circular)?  polariitd.  Tit  incident  plan*  wav*  i*  circularly  polariud. 

Th*  angle*  of  incidence  (t*ta»  ,  ia  degree*,  Tb*  angle*  of  incidence  (t*ta,p)  ,  ia  degma, 

labeling  tb*  r„  curvet  an  related  to  lb*  con*-  labeling  tb*  r„  curvet  art  related  to  th*  com- 

moa  erientatioa  of  tbt  clutter*.  moa  orientation  of  tb*  clutter*. 


On  the  basis  of  the  results.  shows,  we 
can  conclude  that  the  aborption  efficiency  of  a 
random  dispersion  of  identical  scattc.crs 
appreciably  depends  om  the  degree  of  anyso- 
tropy  of  tha  single  scatterer*.  Within  large 
limits  of  the  site  parameter,  however,  several 
kinds  of  aaysotropy  give  very  similar  results, 
as  far  as  the  compactness  is  comparable.  Ulti¬ 
mately,  in  a  random  dispersion  the  site  distri¬ 
bution  of  the  ecatterert  aa  inferred  from  opti¬ 
cal  data  should  not  be  considered  reliable  if 
tbe  single  scattcrers  are  appreciably  aaisotro* 
piel**. 
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Density  dependence  of  the  absorption  coefficient 
of  a  dispersion  of  spherical  metal  particles* 

E.  Fucile,  P.  Denti,  R.  Saija,  F.  Borghese 
Universita  di  Messina,  Istituto  di  Struttura  della  Materia 
Salita  Sperone  31, 98166  S.  Agata-Messina,  ITALY 

In  this  paper  we  make  an  attempt  to  assess  to  what  extent  the  absorption  coefficient 
of  a  dispersion  of  spherical  scatterers  is  influenced  by  the  multiple  scattering  effects.  To 
this  end  the  scatterers  of  the  dispersion  are  considered  as  aggregated  in  pairs  (even  if  they 
are  not  in  contact).  The  field  scattered  by  each  pair  is  the  computed  by  means  of  a  the¬ 
ory  we  developed  in  the  last  few  years  just  to  deal  with  aggregates  of  spherical  scatterers. 
Of  course,  the  field  scattered  by  each  pair  depends  on  its  orientation  with  respect  to  the 
incoming  field  as  well  as  on  the  separation  of  the  component  spheres.  By  means  of  the 
theory  cited  above,  we  are  able  to  sum  analitically  over  the  orientations  to  get  the  a  scat¬ 
tered  field  that  is  still  a  function  of  the  distance  of  the  spheres  composing  the  pair.  At  this 
stage  we  make  an  weighted  average  over  the  distances  using  as  a  weight  the  pair  correla¬ 
tion  function  for  a  dispersion  with  the  density  we  are  dealing  with.  As  a  result  we  get  the 
field  scattered  by  each  sphere  dressed  by  the  effect  of  multiple  scattering  processes  pro¬ 
duced  by  the  other  particles.  Of  course  this  procedure  yields  only  a  first  approximation  as 
only  pair  effects  are  included  while  those  produced  by  triplets  and  higher  order  multiplets 
of  particles  are  neglected.  We  are  able  to  show  that  the  results  we  get  amount  to  solve 
iteratively  the  exact  Foldy-Twersky  integral  equation  up  to  the  third  order  iteration. 


•Based  on  work  supported  in  part  by  the  U.  S.  Army  European  Research  Office  through  Contract 
DAJA45-86-C-0003  and  in  part  by  the  Consiglio  Nazionale  delle  Ricerchc  through  the  Gruppo  Nazionalc 
Struttura  della  Materia 
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It  is  well  known  that  the  most  common  aerosols  are  composed  either  of, 
irregularly-shaped  particles,  or  of  spherical  particles  that  may  undergo 
aggregation  phenomena  which  in  their  turn  give  rise  to  effectively  non  iso¬ 
tropic  compound  scatterers.  It  is  therefore  not  surprising  that  the  inter¬ 
pretation  of  the  optical  behaviour  of  aerosols  based  on  the  well  known  Mie 
theory  are  often  unsatisfactory. 

A  few  years  ago  we  devised  a  method  suitable  to  calculate  the  scatte¬ 
ring  of  electromagnetic  waves  form  a  cluster  of  spherical  scatterers  of 
known  geometry.  Our  purpose  was  to  account  for  the  effect  of  the  anysotropy 
of  the  scattering  particles  without  undue  computational  effort  and  to  see 
to  what  extent  a  truely  irregularly-shaped  particle  could  be  approximated 
by  a  simpler  system. 

To  describe  the  scattering  by  such  a  compound  scatterer  we  Hiake  full 
use  of  the  possibility  of  expanding  both  the  incident  and  the  scattered 
field  in  a  series  of  spherical  multipoles.  As  a  consequence  we  are  able  to 
describe  both  the  electric  and  the  magnetic  multipoles  on  the  same  footing 
without  introducing  any  approximation  but  the  truncation  of  the  multipolar 
expansions  in  order  to  get  a  finite  set  of  equations. ( 1  ) . 

The  calculations  we  performed  on  clusters  of  various  geometry  show, 
as  expected,  a  noticeable  dependence  on  the  extinction  cross  section  (as 
well  as  of  the  scattering  cross  section)  on  the  polarization  and  on  the 
direction  of  incidence  of  the  plane  wave.  Furthermore,  thank  to  the  trans¬ 
formation  properties  of  the  spherical  multipoles  under  rotation,  we  were 
able  to  show  that  the  forward-scattering  amplitude  of  any  cluster  depends 
analytically  on  the  direction  of  incidence  in  a  rather  simple  way. 

All  the  features  mentioned  above  are  of  invaluable  help  when  comparing 
the  results  of  our  approach  to  the  experimental  data.  Nevertheless,  before 
we  make  the  actual  comparison  with  the  experimental  measurements  of  Schuer- 
manv  Wang  and  Greonberg( 2  )  we  recall  some  facts  about  the  convergency  of 
any  calculation  based  on  multipolar  expansions. 

For  a  single  sphere,  of  radius  b,  the  convergency  of  the  calculations 
depends  on  the  product  x=kb,  the  so  called  size-parameter  of  the  scatterer, 
k  being  the  magnitude  of  the  incident  wave  vector.  If  L^  is  the  maximum 
multipolar  order  entering  the  calculations,  the  convergency  requires  LM>x. 
For  a  cluster  of  spheres  we  showed  that  the  parameter  corresponding  to  x  is 
xo=kbo  •  bo  bein9  the  radius  of  the  smaller  sphere  that  can  include  the 
whole  cluster.  Even  in  this  case  a  good  convergency  is  achieved  only  for 
lM>x0.  In  any  case  scattering  objects  with  the  same  x  (or  x„l  have  the  same 
scattering  features,  provided,  of  course,  that  the  refractive  index  is  non 
dispersive . 

In  figure  1  we  show  the  experimental  setup  used  by  Schuerman  and  Wang 
to  perform  their  mesurements  on  single  clusters  as  a  function  of  the  direc¬ 
tion  of  incidence.  We  notice  that  to  ensure  the  reproducibility  of  the 
results  only  the  polarization  indicated  in  figure  1  can  be  used. 


S phf  A 


Tn  figure  2  and  3  wo  report  the  results  of  our  calculations,  for  both 
orthogonal  and  parallel  polarization,  for  a  cluster  of  two  spheres  as  a  fun 
ction  of  the  direction  of  incidence,  together  with  the  corresponding  experi 
mental  results  of  Schucrman,  Wang  and  Greenbergi  3  ) .  More  precisely,  on  the 
axes  are  reported  the  real,  0,  and  the  immaginary,  P,  parts  of  the  forward 
scattering  amplitude  of  the  cluster. 


Fig.  2 


Fig.  3 


o  Orthogonal  polarization  (aiparlmantal,  ral.3) 
*  Orthogonal  polarization  (calculator!) 

□  Paratlal  polarization  (calculator!) 


Both  figure  2  and  figure  3  show  that  the  agreement  between  the  experimental 
and  the  theoretical  results  is  rather  good.  We  state  this  also  on  the  folio 
wing  considerations.  First,  Schuerman  and  Wang( 3 )  extimate  the  experimental 
error  in  about  KMC  in  magnitude  and  12®  in  phase.  Second,  the  angular  depen 
donee  of  the  forward  scattering  amplitude  can  be  determined  analytically, 
and  our  theoretical  curves  are  consistent  with  such  angular  dependence. 

We  want  only  to  stress  that  according  to  the  parameters  reported  in 
fig. 2  and  fig. 3  the  convergency  of  our  calculations  would  require  at  least 
Lm=6.  Actually  we  have  to  use  LM=10  for  both  cases.  Any  further  increase  of 
Lk  has  no  effect  on  the  results  thus  proving  that  we  get  the  best  converge- 
ce  achievable  through  this  kind  of  expansion  of  the  field. 


Fig.  4 


Head-on  Incidence  (numbers  give  ks) 
o  experimental,  ret.  2 


•  calculated 

■  (Q,P)  as  calculated  tor  two 
Independent  spheres 


In  figure  4  we  report  the  results  for  a  cluster  of  two  spheres  whose 
distance  is  gradually  increased.  On  account  of  the  considerations  reported 
above,  when  x0  the  convergency  would  require  Ljn?30.  This  cannot  be  done 
without  using  delicate  numerical  procedures  which  in  turn  require  a  very 
large  and  fast  computer.  For  this  reason  our  present  results  fit  satisfac¬ 
torily  the  experimental  results  only  for  small  values  of  x0,  although  the 
uncertainty  of  the  experimental  data  must,  even  in  this  case,  be  taken  in¬ 
to  account. 

In  conclusion  we  can  state  that  our  approach  is  able  to  describe  cor¬ 
rectly  the  scattering  from  a  cluster  of  spheres.  At  present  we  are  impro¬ 
ving  the  programs  to  include  the  possibility  of  extending  the  multipolar 
expansion  to  very  high  values  of  L^.  In  any  case  both  the  theoretical  and 
the  experimental  approach  mutually  support  the  respective  findings. 
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